More efficient Bell inequalities for Werner states 



00 

o 
o 

(N 
(N 



Oh' 
I ■ 

> • 

■ 

a ■ 

=3 ; 



> 

o\ 
o 
o 

\6 
o 
oo 

o 



T. VcrtcsE 

Institute of Nuclear Research of the Hungarian Academy of Sciences 
H-4001 Debrecen, P.O. Box 51, Hungary 
(Dated: June 2, 2008) 

In this paper we study the nonlocal properties of two-qubit Werner states parameterized by the 
visibility parameter < p < 1. New family of Bell inequalities are constructed which prove the two- 
qubit Werner states to be nonlocal for the parameter range 0.7056 < p < 1. This is slightly wider 
than the range 0.7071 < p < 1, corresponding to the violation of the Clauser-Horne-Shimony-Holt 
(CHSIf) inequality. This answers a question posed by Gisin in the positive, i.e., there exist Bell 
inequalities which are more efficient than the CHSH inequality in the sense that they are violated 
by a wider range of two-qubit Werner states. 

PACS numbers: 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

Quantum Mechanics is inherently nonlocal, clearly 
demonstrated by the fact that measurements on quantum 
states may violate the so-called Bell inequalities [H, S]. 
This has been verified experimentally as well, up to some 
technical loopholes On the other hand, when a quan- 
tum state cannot be prepared using only local operations 
and classical communication, it possesses quantum cor- 
relations and we say that the state is entangled. It was 
Werner who asked firstly what the relation is between 
quantum nonlocality and quantum correlations It is 
actually known that any pure entangled state of two or 
more subsystems may violate a generalized Bell inequal- 
ity 0, Q , thus here nonlocality and entanglement coin- 
cide. For mixed states, however, the relation between 
entanglement and nonlocality is much complicated. In 
1989 Werner Q constructed a family of bipartite mixed 
states (became known as Werner states), which, while be- 
ing entangled, yield outcomes that admit a local hidden 
variables (LHV) model. This conclusively proved that 
entanglement and nonlocality are different resources. 

However, if we want to describe quantitatively the dif- 
ference, the picture turns out to be quite subtle even in 
the case of two-qubit Werner states, which are mixtures 
of the singlet = (|01) - |10))/V2 with white noise 
of the form 



p7 



p|^-)(7^-| + (l-p)l/4. 



(1) 



Werner showed [3l that these states are separable if and 
only if p < 1/3. With respect to the locality properties, 
on one hand, Werner states admit a LHV model for all 
measurements for p < 5/12 Q and admit a LHV model 
for projective measurements for p < 0.6595 Q. On the 
other hand; Werner states violate the CHSH inequality 
for p > l/\/2, in which case LHV model clearly cannot 
be constructed. It is not known whether Werner states 
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admit an LHV model in the region 0.6595 <p< I/V2. 
The actual value of p where the state ceases to be non- 
local, designated p^ , is particularly relevant from the 
viewpoint of experiments since this value specifies the 
amount of noise the singlet tolerates before losing its non- 
local properties. This issue was addressed by Gisin some 
time ago @ (see also [l3|), who posed the question to 
find Bell inequalities which are more efficient than the 
CHSH one for Werner states. In this paper we intend to 
give a definite answer to this question by providing Bell 
inequalities which can be violated slightly stronger than 
the CHSH one, resulting in the bound p^ < 0.7056 for 
the nonlocality visibility threshold (instead of the bound 
< I/V2 ~ 0.7071 owing to the CHSH inequahty). 
Note, how powerful the CHSH inequality is, which is the 
simplest Bell inequality, consisting of two settings on each 
side, while a stronger Bell inequality presented in this 
work has at least 465 settings on each side. 

We also would like to point out that while in certain 
cases using sequence of measurements may extend the 
range of locality [ll|, the nonlocality threshold p^ for 
Werner states could not be decreased even on this way 
[1^ . There is also an interesting line of research, which 
explores the parameter region of Bell violation for Werner 
states by restricting the class of possible LHV models 

El, 

Actually, Ref. [l3| could achieve violation of 
certain Bell inequalities, assuming the above limitations 
for p > 1/3, i.e., for the entire range of the nonsepara- 
bility region. An other way of generalization to obtain 
the range of locality is the extension of the Werner states 
to e.g., more parties 15 1 or higher dimensions p^. [l7|. 
However, let us mention, that a gap also remained in 
these cases between the best known local model [l^ [l3| 
and the proven nonlocality threshold [11, [III ■ 

The outline of the present work is as follows. In Sec.HIl 
we briefly summarize the relation between Bell inequal- 
ities for two-qubit Werner states and Grothendieck con- 
stant of order 3, denoted by KciS). In Sec. IIIII a family 
of Bell inequalities is constructed and in Sec. IIVI with the 
aid of these inequalities a lower bound, bigger than v^, 
is given for Kc{5), implying that Werner states ([T|) with 
p < l/\/2 can still violate these inequalities. In Sec. |V] 
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we give a better lower bound for iff? (3) and for higher 
orders {Kcid) with d = 4, 5), as well In Sec. |Vl]a Bell 
inequality is provided with a number of settings 11 and 
14, proving Kg{4:) > v^, and in Sec. IVIll the relevance 
property of the constructed family of Bell inequalities 
is demonstrated. Sec. IVIIII summarizes and poses some 
open questions. 

II. BELL INEQUALITIES LINKED TO 
GROTHENDIECK CONSTANTS 

Define the expression 

m 

I ^ \ J2 a.b,\, (2) 

where M is any m x m matrix with real entries and 
ai, . . . , flm, 6i, . . . , bjn G {^1, +!}■ Now let us define 

m 

/•^^l ^ M,,a,.6,|, (3) 

where the unit vectors oi, . . . , dm, bi, . . . , bm are in M'' 
and a • & is the dot product of a and b. Grothendieck 
constant plays a prominent role in the theorem of linear 
operators on Banach spaces [131 ■ Grothendieck constant 
of order d, designated Koid), for any integer d > 2, can 
be defined as [21| 

/'^ < Raid) max / (4) 

ai,bj—±l 

for all unit vectors ai , . . . , , fei , . . . , in K'' and for 
all m X m matrix M. The constant Koid) is taken to be 
the smallest possible one. 

Now let us discuss briefiy the connection with Bell in- 
equalities. In the Bell scenario we consider two parties, 
Alice and Bob, each chooses from m ±1- valued observ- 
ables, specified by {Ai, . . . , Am} and {Bi, . . . , Bm}- The 
joint correlation of Alice and Bob's measurement out- 
comes, designated ai and (3j respectively, is given by 
{ai(3j) = Tr{Ai (8) Bjp), where p denotes the density ma- 
trix of the bipartite state. A correlation Bell inequality 
can be written as 

m 

M,,{a,(3,)<L, (5) 

where L signifies the bound which can be achieved by 
local models and M is a m x m matrix with real coef- 
ficients defining a Bell inequality. The local bound can 
always be achieved by a deterministic local model, i.e., 
for all real numbers a^, 6j = ±1 we have 

m 

max y MijGibj = L. (6) 
IJ=1 



In this way the expression / defined by ([2]) is linked to 
a correlation Bell inequality with matrix M and local 
bound maxa- f,j=±i I ~ L. 

On the other hand, for the singlet state p = 
we have {ai/3j)^- = (8) Bj\il;~) = —di ■ bj, where 

the observables A — dd and B — ba corresponding to 
Alice and Bob's projective measurements are specified 
by the unit vectors a and b in K^. Then substitut- 
ing into ([5]) one obtains the expression in ([3]). Fur- 
thermore, Tsirelson [2^ proved that correlations which 
are dot products of unit vectors d,b £ M.'^ can always 
be realized by performing projective measurements on 
maximally entangled states in some higher dimensional 
Hilbert spaces. Thus the value max/'' can always be 
achieved by means of quantum mechanics. 

Since joint correlations vanish for the maximally mixed 
state, it follows that the critical point at which Werner 
states in (jlj cease to violate any Bell inequality is = 
1/A'g'(3). This key correspondence has been established 
in Ref. [Sj]. Though, the exact value of A'g'(3) is not 
known, but known bounds establish that 0.6595 < < 
0.7071. In this paper we show that Ag'(3) > 1.4172 im- 
plying the shghtly smaller gap 0.6595 < p^ < 0.7056. 
Let us mention, that the Fishburn-Reeds Bell inequality 
[2^ provides an explicit example with 20 settings on each 
side showing that A'g(5) > 10/7 = 1.42857 > V2. Also 
Toner has shown that Ag(4) > [§]. But, as far as we 
know, the question has remained open whether Kg{3) is 
bigger than ^2 implying p^ < 1/ V2- 



III. CONSTRUCTING FAMILY OF BELL 
INEQUALITIES 

For Bell diagonal states, such as for Werner states, 
under projective measurements Alice and Bob's local 
marginals (defined by (ai) = Tr{Ai (g) Ip) for Alice and 
likewise for Bob) are zero, thus it is sufficient to deal 
with generic correlation Bell inequalities defined by ^ to 
obtain maximal Bell violation for Werner states. More- 
over, in this respect, the tight correlation Bell inequali- 
ties, which can be considered as facets of the correlation 
polytope [i^l, specified by the number of two-outcome 
measurements m on each side, are the most efficient ones. 
For TO = 2 one obtains as the only nontrivial correlation 
inequality the CHSH one [2^. For to > 2 one needs to 
resort to numerical programs for computing the inequali- 
ties corresponding to the inequivalent facets of the corre- 
lation polytope. Up to TO = 4 all the correlation inequal- 
ities have been computed [26j . and the two inequivalent 
inequalities obtained are in fact less efficient than the 
CHSH one for Werner states. However, the complexity 
of the computation exponentially grows with to (in fact, 
this is an NP-complete problem [2J|), therefore there is 
no hope to completely characterize all the facets of the 
correlation polytope for any given to. Thus in general one 
needs to look for alternative methods. For instance Gisin 
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explored special form of families of tight correlation in- 
equalities, the so-called D-inequalities in Ref. Avis 
et al. [2^1 applied triangular elimination to the list of 
known facet inequalities of the cut polytope to construct 
many new tight correlation inequalities. Alternatively, 
one can construct (possibly not tight) correlation inequal- 
ities which however can be easily generalized to arbitrary 
number of settings, such as in the cases [13, |28,, 29,, 30| . 
In the present work we have chosen this latter direction 
by modifying the correlation inequalities Z„ introduced 
in H^. 

Let us specify the form of M in ([2]) through the follow- 
ing formula. 



i=l j = l 

l<.i<ij<.nB l^i<J^''T'A 

(7) 

entailing a Bell inequality with niA = UA + nsinB — l)/2 
and niB = TLB -\- nAin-A ~ l)/2 measurement settings on 
Alice and Bob's respective side. First we calculate the 
maximum achievable value (local bound) for it. For this 
sake, we can write for the maximum 



max In A ,„B = max { Oi bj 

ai ,bi—±l ' ^ — 



l<i<j<nB 



l<'t<j<7l^ 



(8) 



Generally, the local bound can be obtained by finding the 
maximum over all possible values ai,bj = ±1. However, 
in this particular case one can exploit the symmetry with 
respect to change of indices within the sets {ai}"^^ and 
Thus one needs to check altogether uaub cases 
where -1-1 occurs 1 < k < ua times in the set {ai}"^^ 
and occurs 1 < I < ub times in the set {b}^^^ (the 
rest being —1). For any k,l pair we have max/„^.„g = 
max{(nA - 2k){nB - 21) + 2{nA - k)k + 2{nB- 1)1}. This 
expression is maximal by fc — / = [(n^ — 77,b)/2J resulting 
in the local bound 



in the two parties, reads 

n n 

itn^YY'^^-^^ 

i 3 

+ Y a,j-{b,-bj)+ Y hj-{ai-aj). 

(10) 

For the maximum, similarly to the LHV case, the two- 
indices terms can be omitted: 

max J^_„ = max{^ a^Y^^J + Y ' ~ ' 

^ n 1 ^ 

+ Y <max{-|^a,p + -|^fe,|2 

l<i<j<n i— 1 i—1 

+ Y Y (11) 



l<i<j<n 



l<z<j<n 



where the maximization is over all a;, bj G S'^^^ and the 
last inequality comes from the relation between the geo- 
metric and quadratic mean. Furthermore, since {fli}"^]^ 
and {biYl^i '^^ ^ot depend on each other, one can maxi- 
mize the two sets independently resulting in 

n 

max/^ „ = max{| Y + 2 Y " "j'^' (1^) 

with the constraints Si £ S"^^^, where the equality sign 
is due to the fact that at the maximum one can take 
di — hi for all i, which saturates the inequality in (jlip . 
This expression shows some similarity with the one ap- 
pearing in [soj . in which case one had to maximize only 
the last term, i.e., the sum of distances of n unit vectors. 
This is a problem occurring in discrete mathematics, and 
there exist optimal solutions for various instances of the 
values n,d [3ll |. In contrast, in the present case, the 
quadratic term makes things complicated to get the true 
optimum value. However, due to formula ^ one can 
give the lower bound KQ{d) > /^^ max/,i^_„g for 
the Grothendieck constant of order d, without knowing 
the true maximum value max /f . „ „ . 



IV. LOWER BOUND FOR GROTHENDIECK 
CONSTANT OF ORDER 3 



max/„^^„g =("i + l)/2, for \nA - ub] odd 
max/„^_„g =("-A + '^b)/2, for \nA - even. (9) 

In this paper we focus on two particular cases ua = 
ns-l-l and ua = riB, but first let us restrict our attention 
to the latter, symmetric case ua — nB — n. The LHV 
bound gives max/„.„ = m? by inserting ua ~ riB = n in 
0. On the other hand, the expression /^„, symmetric 



In fact, in the particular case /"„ one can obtain the 
exact maximum, max/"„ = 3/2 — 1/(272). This result 
(noticing that for large n the violation tends to 1.5) would 
indicate that there may be some hope to get a lower 
bound Kg{^) > \f2. Below we show that the maximum 
above can indeed be attained. 

First let us observe that in (|12p only n vectors occur, 
thus we have majil^n = max/^„ with m = n{n + l)/2, 
where m is the number of measurement settings on each 
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side. Now we take the unit vectors Ei in such a way 
that Hi ■ dj = 1/2 for all i ^ j- This can be achieved, by 
noting that the nxn Gram matrix G, defined by elements 
Gij = di-dj is positive definite, and every positive definite 
matrix is a Gram matrix for some set of vectors. Thus 
it is enough to show that the Gram matrix G, defined 
as above {Gij = 1, Vi = j and Gij = 1/2, Vi ^ j), is 
positive definite. However, using the Sylvester criterion 
[3^ one can establish that G defined above is positive 
definite iff det G > for any n. One may obtain by 
induction the closed formula det A = (a — 6)"^^ (a + (n — 
1)6) for the determinant of any nxn matrix A having 
in the diagonals the value a and in all off-diagonals the 
value b. By choosing particularly a = 1 and 6=1/2 one 
gets det G > for any dimension n, which proves our 
assertion. 



Beside, all the elements in (|12p can be obtained as the 
only function of di ■ dj, since | Yl7=i ~ '^i ' % 

and \di — dj\ = yj2 — 2di ■ dj . Thus by substitution we 
obtain „ = n{3n — 1)/2. Then it follows using ^ , that 
/"„/max/„_„ = 3/2 — l/(27i). It is also possible to verify 
that this is in fact the maximum value. The verification, 
which is not detailed here, goes the same line as discussed 
by Wehner in [s^ for the chained Bell inequality p7| 
through the dual solution of a semi-definite optimization 
problem [s^ . Note that this optimization problem is just 
the first step in the hierarchy introduced by Navascues 
et al. fssl. [36j. 

Now we wish to obtain a lower bound Kg(3) > 
„/ max/„^„ bigger than \/2 for the Grothendieck con- 
stant of order 3, Kg{3)- Since owing to (O, max/„.„ = 
n^, we are left with the calculation of /^^ „ which though 
might be not maximal, but large enough to supply a good 
lower bound for Kg{3). This is achieved by substituting 
in in the place of di explicit values on the follow- 
ing way. Since di are unit vectors in M."^, the first two 
components Pi = {xi,yi) of di, which can be considered 
as points in the XY-planc, completely specify the vec- 
tor itself. Let n = 30 and distribute these points on 
three co-centric circles centered at the origin (0, 0) with 
radii pi = 22/100, pu = 52/100 and pm = 77/100. 
Then let Pi — (p/ cos 7r/4, p/ sin7r/4), P5 = (0,^//) and 
P15 = {0,piii)- The other Pi vectors arc constructed 
from the above vectors by rotating them with angles 
7r/2, 7r/5, and 7r/8, respectively, such as to form regu- 
lar polygons with vertices 4, 10 and 16 (as it is shown 
in the figure). By inserting the explicit values of the 
corresponding set {di}i2^i into the expression to be max- 
imized in (fT2l) one obtains I^„/In,n — 1.415 199 with 
n = 30. This implies that Kg{3) is indeed bigger than 
y/2 = 1.414 2136 . . .. The specific values of di have been 
found by performing optimization with respect to the 
radii of the three circles by choosing regular polygons 
with various number of vertices. 
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FIG. 1: The 30 points which are projection of the vectors 
Si on the XY -plane. They are equally distributed on three 
co-centric circles with radii 22/100,52/100,77/100 centered 
at the origin. The outer circle represents the grand circle 
projected on the XY-plane, thus having radius 1. 



V. BETTER LOWER BOUNDS FOR 
GROTHENDIECK CONSTANT OF ORDERS 3,4,5 

In this section a general method is discussed to obtain 
local maximum on if^ „ for any n, d, which for many in- 
stances are presumably the global or close to the global 
maximum. Then, recalling KQ^d) > /^„/max/„.„ and 
max/„_„ = n^, this method yields lower bounds for 
Kcid). In particular we present results for d = 3,4,5, 
calculated by the value n ~ 100. 

Let us consider the following iteration scheme, which 
is a simplified version of the see-saw iteration method, 
already used in the literature to solve optimization prob- 
lems in similar context entering many optimization pa- 
rameters [l^jll^ . Note, that the matrix M of defined 
through ([7]) is symmetric, thus we may write 

m mm mm 

(13) 

with m = 7i(?i + l)/2 and di,hj are unit vectors in 
. In this notation we contracted the double indices 
ij appearing in pU]) . so that {di}^]^ stands for the set 

({^i}"=ii {aij}i<i<j<n) and similarly for the vectors h. 
Considering p3p one can maximize the expression 
„ for given {a;}™]^ by setting hi parallel to 
for all i. Then one can continue with setting di parallel 
to J2 Mijbj for all i. However, due to the fact that M is 
symmetric, one can get rid of the vectors b and obtain 
the iteration rule di J2j ^-^ij^j /\J2j -^ij'^jl ^O'' 
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provided | MijSjl ^ 0. Here the notation |?7| refers to 
the Euchdean norm of a vector v GMf^. Thus our task is 
to give initial values for the unit vectors di G K'*, which 
we choose in the following way. The surface of the unit 
sphere in Mf^ can be parameterized by d — 1 angles, Si = 
(cos((Ai),sin(</.i)cos(0f), . .^,sin(</>i) • ■ ■ sinicj^t^) cos{q>t') 
sm{4>l) ■ ■ ■ sm{(j)f ^)sin(0f ^)). We define the starting 
vectors {a^}™ ^ with angles 0*'' = k [rad], 1 < A: < d — 1. 
Then we perform the above iteration scheme for a given 
time, practically till the vectors di in two successive 
iteration steps differ less than an infinitesimal thresh- 
old value. In particular, for n = 100 we found that 
1000 iteration steps were sufficient for our purposes. 
Also, we checked for each case d = 3,4,5, that the 
value {"^j Mijdjl in the denominator of the iterated 
expression was nonzero (actually, it was no less than 
10^** for all i in each case of d). On the other hand, 
the iteration was performed with machine precision 
~ 10~^^ in the Mathematica package, and we checked 
that \d, • a, - 1| < 10"^^ for alll < i < m after the 1000 
iteration steps completed. 

For n = 100, we obtained the following numbers, 
/3 „/max/„,„ = 1.417 241, /4 „/max/„,„ = 1.445 207, 
and J^„/max/„_„ = 1.460 065. These numbers are 
lower bounds for the Grothendieck constants Kci'i), 
Kc{4:) and Kg{5), respectively. We mention that for 
n = 100 the dimension of the respective matrix M is 
n(n + l)/2 = 5050. Note that the best lower bound for 
Kg{5) presented so far in the literature Kg{5) > 10/7 = 
1.428 571 . . . comes from the Fishburn-Reeds inequality 
[2^. Our result for /vg(3) provides us with the better 
lower bound pj^ < 0.705 596 for the critical value 
owing to the formula = 1/A'g(3). 



VI. MINIMAL NUMBER OF MEASUREMENTS 

One may also ask, what is the smallest number of set- 
tings on Alice and Bob's side, where Kcid) can exceed 
V2 for some d > 2. To the best of our knowledge, so 
far it has been provided by the Fishburn-Reeds inequal- 
ities [231 . Their construction, giving Kq{5) > 10/7 = 
1.428 571 . . . can be obtained by 20 measurement set- 
tings on each side. Now we choose n^=:nB + l = 5in 
expression InA,nB of giving the number of settings 
11 and 14 on Alice and Bob's side, respectively. Thus 
the matrix M in this particular instance has dimensions 
11 X 14. We show that this expression I^a provides us 
with an example where Kq{4) > \/2. 

Substituting values ha — 5 and ub = 4 into the for- 
mula Q for odd \nA — "si one obtains the value 20 for 
the local bound. On the other hand, the maximum value 
corresponding to the vectorial case max/|4, can be ob- 
tained by the mean of semidefinite techniques [1^ as a 
first step of the hierarchy in [3^, where we used the Se- 
DuMi package [s^ for Matlab by the explicit numerical 
computation. This algorithm solves both the primal and 



the dual optimization problem at the same time and thus 
yields bounds on the accuracy of the obtained solution 
as well. Actually, we obtained the same optimal value 
28.390 139 for both cases. This yields the ratio 1.419 507 
for the violation of the Bell inequality J5.4 < 20 clearly 
^ beating the V2 limit with 11 and 14 settings on Alice 
and Bob's side, respectively. 

VII. TIGHTNESS AND RELEVANCE OF BELL 
INEQUALITIES 

It would be interesting to know whether the family 
of correlation inequalities defined by ([7]) is tight, i.e., 
whether it is a facet or not of the local Bell polytope [40l | 
consisting of local marginals as well. This can be done 
by computing the dimension of the subspace spanned by 
all deterministic strategies saturating the inequality. If 
this subspace is found to be a hyperplane with dimen- 
sion d = mAinB + mA + tub^ then the inequality is 
tight. Numerically, we treated the ua = ub + I and 
"A = i^B = " cases in the expression InA^ns- Computa- 
tionally we found that in the former case the inequality 
is tight up to UA ~ 4. On the other hand, the latter 
symmetric inequality proved to be not tight, but by the 
addition of some local terms ai^hj = ±1 as follows, 

n 

n n 
l<i<^'<n i—1 J — 1 

we checked computationally its tightness up to n = 4. 
Note, that here the terms a;, hj refer to Alice and Bob's 
local marginals in the corresponding Bell inequalities. 

Recalling from Sec. IIIIl that for = we have 
maxJ„^„ = max{7i^ — 2(/s — ^)^}, by adding the marginals 
max/;'„ = maxjn^ - 2{k - If + {2k - n) - {21 - n)} = 
max{n^ — 2{k — l){k — I + 1)} = n^, thus the local bound 
does not change. Let us notice, that 2 specifie d by 
n = 2 in HH) is just the 13322 Bell inequality [13, lilf . 
which is known to be tight [i^]. In both cases, /„.„_i 
and „, we suspect that these Bell inequalities are tight 
for any higher values of n > 4, as well. 

Let us discuss the concept of relevant Bell inequali- 
ties, whose definition we quote from [l3|, Sec. A.l: "An 
inequality is relevant with respect to a given set of in- 
equalities if there is a quantum state violating it, but not 
violating any of the inequalities in the set." Collins and 
Gisin [4l| showed that the 13322 inequality is relevant to 
the famous CHSH inequality Interestingly, they also 
found that given 13322 the CHSH inequality is no longer 
relevant. Furthermore, Ito, Imai and Avis (26j have re- 
cently conjectured supported by numerical optimization, 
that there exist Bell inequalities relevant for the 73322 in- 
equality for 3-level isotropic states. However, limiting the 
Hilbert space dimension to a qubit pair, they did not find 
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a Bell inequality which would be relevant with respect to 
^3322- Our new inequalities, In^n and „ with n = 100, 
however, are examples to this latter case, demonstrating 
that in the parameter range 0.705 596 < p < 0.7071 107 
two-qubit Werner states do not violate the 13322 inequal- 
ities but violate /„.n or 7,'^ „ for n = 100 (note that for 
the Werner states the local marginals become identically 
zero, thus in this respect /„^„ and „ are equivalent). 

Moreover, one can demonstrate, that there is an inclu- 
sion relation, a notion introduced in [4^, between „ 
and In_i meaning that one can obtain the inequality 
In-i n-i measuring the identity for some settings in 
the inequality „ (i.e., performing degenerate measure- 
ments) . This implies that „ for any n > 2 is relevant 
with respect to /2 2 = ^3322. The proof is simple, actu- 
ally by setting a„ = -t-1, 6n = +1 and ai„ = — 1, bin = ^1 
for 1 < i < n in n oi^c obtains /^-i n-ii and then by 
induction one arrives at /2 2 • 

Altogether, one can say that if one limits the Hilbert 
space dimension to two qubits ([l3|, Sec. A. 2) the 1^ ^ 
inequality for n — > cx) is the only relevant one with respect 
to all presently known Bell inequalities. 

VIII. SUMMARY 

We provided a new family of Bell inequalities which 
proves that Werner states in ([1]) arc nonlocal for the 
parameter range p > 0.7056, the best earlier result 



p > 0.7071 is given by the CHSH inequalities. Some 
of these Bell inequalities are shown to be relevant with 
respect to any other known Bell inequality. Our results 
have been obtained by proving that the Grothendieck 
constant of order 3, Kg{3), is bigger than V2, in partic- 
ular, 7^0(3) > 1.4172. Though our result for the wider 
visibility range of nonlocal Werner states has been ob- 
tained for a number of settings (at least 465 for each 
party) which are not particularly suited for experiments, 
we believe that they are interesting from a conceptual 
point of view. Entangled states in many quantum in- 
formation protocols (for instance in quantum communi- 
cation complexity [46| and device-independent quantum 
key distribution problems [i^l) give advantage over there 
classical counterparts only if they exhibit nonlocal cor- 
relations. Thus, in this paper we have shown that this 
nonlocal correlation can in principle be exploited in a 
wider range of Werner states. 

We leave it as an open question how to construct even 
better inequalities which would allow to beat the limit 
of Kci'i) stronger. The possibility for such inequali- 
ties is suggested by the fact that an upper bound for 
Kg = lim„^oo Kcin) is 1.7822 which is suspected to be 
tight [21I, But the inequality for n ^ cxd gives 
the lower bound 1.5 for Kg, which is even smaller than 
the lower bound 1.6770 for Kg presented in Ref. [isj . 
Thus it is not impossible that there exist inequalities pro- 
viding bigger values for 7^0(3) entailing even better Bell 
inequalities than the present ones for Werner states. 



J.S. Bell, Physics 1, 195 (1964). 

J. Clauser, M. Horne, A. Shimony, and R. Holt, Phys. 

Rev. Lett. 23, 880 (1969). 

A. Aspect, Nature 398, 189 (1999). 

R.F. Werner, Phys. Rev. A 40, 4277 (1989). 

N. Gisin and A. Peres, Phys. Lett. A 162, 15 (1992). 

S. Popescu and D. Rohrlich, Phys. Lett. A 166, 293 

(1992). 

J. Barrett, Phys. Rev. A 65, 042302 (2002). 

A. Acin, N. Gisin, and B. Toner, Phys. Rev. A, 73, 

062105 (2006). 

N. Gisin, problem 19 (2003), presented in the web page 

http : //www. imaph. tu-bs .de/qi/problems. 

N. Gisin, arXiv:quant-ph/0702021v2 (2007). 

S. Popescu, Phys. Rev. Lett. 74, 2619 (1995). 

N. Gisin, Phys. Lett. A 210, 151 (1996). 

1. Pitowsky, J. Math. Phys. 49, 012101 (2008). 

C. F. Wildfeuer and J. P. Bowling, arXiv:0708.1973v3 
(2007). 

G. Toth and A. Acin, Phys. Rev. A 74, 030306(R) (2006). 
M.L. Almeida, S. Pironio, J. Barrett, G. Toth, and A. 
Acin, Phys. Rev. Lett. 99, 040403 (2007). 

H. M. Wiseman, S.J. Jones, and A.C. Doherty, Phys. Rev. 
Lett. 98, 140402 (2007). 

D. Collins, N. Gisin, N. Linden, S. Massar, and S. 
Popescu, Phys. Rev. Lett. 88, 040404 (2002). 

D. Kaszlikowski, P. Gnacinski, M. Zukowski, W. Mik- 
laszewski, and A. Zeilinger, Phys. Rev. Lett. 85, 4418 



(2000). 

G. Pisier, Factorization of Linear Operators and Geome- 
try of Banach Spaces, (American Mathematical Society, 
Providence, RI, 1986). 

S.R. Finch, Mathematical Constants, (Cambridge, Cam- 
bridge University Press, 2003), pp. 235-237. 
B.S. Tsirelson, J. Soviet. Math., 36, 557, (1987). 
P.C. Fishburn and J.A. Reeds, SIAM J. Discrete Math. 
7, 48 (1994). 

I. Pitowsky, Quantum Probability, Quantum Logic, 
(Springer- Verlag, 1989). 

B.S. Tsirelson, Lett. Math. Phys., 4, 93 (1980). 
D. Avis, H. Imai, T. Ito, J. Phys. A 39, 11283 (2006). 
S. Braunstein and C. Caves, Ann. Phys. (NY), 202, 22 
(1990). 

N. Gisin, Phys. Lett. A, 260, 1 (1999). 

H. Bechmann-Pasquinucci and N. Gisin, Phys. Rev. A 
67, 062310 (2003). 

T. Vertesi and K.F. Pal, Phys. Rev. A 77, 042106 (2008). 
H. Cohn and A. Kumar, J. Am. Math. Soc, 20, 99 
(2007). 

R.A. Horn and C.R. Johnson, Topics in Matrix Analysis 
(Cambridge University Press, Cambridge UK, 1991) 
S. Wehner, Phys. Rev. A, 73, 022110 (2006). 
S. Boyd and L. Vandenberghe, Convex Optimization, 
(Cambridge University Press, 2004). 
[35] M. Navascues, S. Pironio, and A. Acin, Phys. Rev. Lett. 
98, 010401 (2007). 



[20 



[21 

[22' 
[23^ 

[24; 

[25' 
[26' 
[27 

[28' 
[29; 

[30 
[31 

[32; 

[33 
[34; 



7 



[36] M. Navascues, S. Pironio, and A. Acin, arXiv:0803.4290 
(2008). 

[37] R. F. Werner and M. M. Wolf, Quantum Information 
& Computation 1, 1 (2001), arXiv:quant-ph/0107093 
(2001). 

[38] T. Ito, H. Imai, D. Avis, Phys. Rev. A 73, 042109 (2006). 
[39] J. Sturm, SeDuMi, a MATLAB toolbox for 

optimization over symmetric cones, URL 

http : / / sedumi .mcmaster . ca. 
[40] M. Froissart, Nuovo Cimento B 64, 241 (1981). 
[41] D. Collins and N. Gisin, J. Phys. A: Math. Gen. 37, 1775 

(2004). 

[42] D. Avis and T. Ito, Discrete Appl. Math., 155, 1689 
(2007), arXiv:quant-ph/0505143 (2005). 



[43] D. Avis, H. Imai, T. Ito, Y. Sasaki, J. Phys. A. 38 10971 

(2005), arXiv:quant-ph/0505060 (2005). 
[44] J. Krivine, Adv. Math. 31, 16 (1979). 
[45] A.M. Davie (unpublished note, 1984) and J. A. Reeds (un- 
published note, 1991). 
[46] C. Brukner, M. Zukowski, and A. Zeilinger, Phys. Rev. 

Lett. 89, 197901 (2002). 
[47] J. Barrett, L. Hardy, and A. Kent, Phys. Rev. Lett. 95, 

010503 (2005); A. Acin, N. Gisin, and L. Masanes, Phys. 

Rev. Lett. 97, 120405 (2006); A. Acin, N. Brunner, N. 

Gisin, S. Massar, S. Pironio, and V. Scarani, Phys. Rev. 

Lett. 98, 230501 (2007). 



